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Abstract 

Strong electron correlations are responsible both for the insulator ground state of undoped 
La2Cu04 and strong antiferromagnetic coupling J between neighbouring spins. We consider mag- 
netic mechanism of superconducting pairing in the effective low energy t — t' — t" — J* model with all 
parameters calculated ah initio. Interaction of strongly correlated electrons with different phonon 
modes is also incorporated. In a BCS type theory the d x 2_ y 2 gap is given by a sum of magnetic and 
phonon contributions. The phonon coupling parameter A = f(x)G, where G is a combination of 
bare electron-phonon couplings for all modes and the function / depends on the hole concentration 
x due to strong electron correlations. The main contribution to the only fitting parameter G is 
determined by a competition of the breathing and buckling modes. Fitting the parameter G from 
the isotope effect we obtain that magnetic and phonon contributions to the critical temperature 
T c work together and are of the same order of magnitude. 
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I. INTRODUCTION 



The discovery of superconductivity in cuprates with high critical temperature T c and 
ci-wave order parameter has serious implications for the theory of superconductivity. The 
presence of large Coulomb interactions that have the potential to destroy conventional s- 
wave BCS state and to transform the half-filled single electron band in undoped La2Cu04 
into the Mott-Hubbard insulator has prompted the search for new approaches for the normal 
and superconducting phases. A microscopic approach based on the conventional ab initio 
local density approximation (LDA) is not valid in the underdoped region of the cuprate 
phase diagram due to the strong electron correlations (SEC) effects In general the 
problem of strong electron correlations is not solved. There are many different theoretical 
approaches to this problem including Quantum Monte Carlo, Exact Diagonalization for finite 
cluster, Slave Boson, Cell Perturbation Theory, Dynamical Mean Field theory for infinite 



lattice 
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The multielectron generalized tight-binding (GTB) method has been proposed 9] which 
describes electrons as the quasiparticle excitations between the local multielectron config- 
urations with the interatomic hopping resulting in the dispersion and the band structure. 
The GTB approach has been successful for undoped cuprates band structure 10J. The 



hybrid LDA+GTB method 



111 ] takes the advantages of both ab initio single electron and 



multielectron approaches. The low energy effective Hamiltonian generated by the quasipar- 
ticle electronic structure is given by the t — t' — t" — J* model (here star means that the 
three-cite correlated hopping ~ J is included) with all parameters calculated ab initio 11]. 
The normal state of the underdoped cuprate is characterized by electron hopping in the spin 
liquid background. At small doping in the antiferromagnetic phase hole is a spin polaron 
and at larger dopin g sp in fluctuations of the short magnetic order also strongly renormalized 



hole dispersion 



n 
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14J. A spin fluctuation pseudogap has been found both in the spin- 



fermion model 15J], in the doped Mott-Hubbard insulator by the cellular dynamical mean 
field theory (DMFT) [l6], and in recent LDA+DMFT+S k calculations 17] for Bi2212. 

A self-consistent consideration of the electronic structure and spin correlation functions 
within the t — t' — t" — J* model with parameters of La2-xSr x Cu04 results in the doping 
evolution of the hole Fermi surface from small pockets around (tt/2, it/ 2) at x < x cr to large 
hole surface around (n, it) at x > x cr with quantum phase transition 18] at x cr ~ 0.15. In 
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this paper we extent this approach to the superconducting phase taking into account both 
magnetic pairing within the t — t'—t" — J* model and electron-phonon interaction (EPI). In 
spite of large number of different phonon modes only a few of them have large EPI including 
the apical oxygen breathing mode (apical oxygen displacement perpendicular the CuC>2 plane 
which modifies the Madelung energy), the in-plane oxygen breathing/half-breathing mode 



19 
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and the buckling mode with in-plane oxygen ions moving perpendicular to the plane 
Recent ab initio study of the electronic structure and EPI in La2- x Sr x Cu04 has proved that 
these three modes contribution to the hole self-energy is more than 80% of the total self- 



energy [21]. 



The multiband p<i-model at low energies is reduced to the effective Hubbard model with 
two Hubbard subbands. The lower and upper Hubbard bands (LHB and UHB) are the 
bands of the Hubbard fermions created by the X-operators Xj a and XJ a ' 2 , correspondingly, 
Xj 9 = \p) (q\. Here |0) , \a = ±1/2) and |2) are the multielectron eigenstates of the CuOq 
unit cell, corresponding to configurations d 10 p 6 with number of holes = 0, d 9 p 6 + d 10 p 5 
with rih = 1, and d 8 p 6 + d 9 p 5 + <i 10 p 4 with nh = 2. In the hole language the electrons in the 
valence band correspond to the holes in the UHB. The amplitudes of interatomic hopping 
in the LHB and UHB are t°® g and t\ g , while the interband hopping is given by tf g . When 
we eliminate the interband excitation throw the charge transfer gap U e /f by the standard 
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unitary transformation the effective Hamiltonian for holes in the UHB is given by 
H t _j* = H t _j + H(_ 3 ) with 

H t -j =E( £ - p)xr + e #*r**r a + 

fcr fgcr 

J2 J fg S g -lrifn g ), 
fg 

fOl +01 
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Where J = ioi/^e// is the super exchange interaction. Distance dependent hopping param- 
eters t^g have been calculated up to 6-th neighbours and it was revealed that contributions 



of the fourth and more distant neighbours are neigligible small 11[. This is the microscopic 
justification of the t — t' — t" — J* model with 3 hopping parameter (t = 0.932, t' = —0.120, 
t" = 0.152, J = 0.298, J' = 0.003, J" = 0.007), all parameters are in eV). The last term H {3) 
corresponds to the three-site correlated hopping that has the same order as the exchange 



term J and has to be included in the theory of superconductivity 



24]. 



In the strong electron correlation regime the EPI is the interaction of phonon (with wave 



number q, frequency ouqy, and mode v) and Hubbard fermions 25|. The effective total 
Hamiltonian is given by 

Heff — H t _j* + H el 

—ph—el (1) 

where effective electron-electron interaction mediated by phonons is H e i_ p h_ e i = 
£ 2 ^kk'q^k+q^k'-q^k' ' 2 X^' 2 with the effective interaction neglecting retardation effects 

kk'q era' 

given by VW q = (K q) 9v ( k '; - q)/^q,w 

For isotropic s-wave gap all phonon modes contributed additively to pairing. For 
anisotropic d x 2_ y 2-g&p the wave number dependence of the EPI matrix elements is crucial. 
Maximal EPI for the breathing/half-breathing mode at large q ~ ir/a results in depairing 
effect of this mode, while the buckling mode with maximum of interaction at q = supports 
the d x 2_ y 2 pairing. This conclusion has been obtained by different approaches 19|, l20l. l26l. |27| 
and results from a simple physics: large q EPI changes the phase of the d^.^-gap on the 
Fermi surface while small q EPI does not change the phase. 

It should be emphasized that mean-field theory in the framework of GTB and GTB+LDA 
methods differs from standard mean field treatment of the Coulomb interaction like Unin 2 —* 
Un-y (77-2), where n, is a number of particles in % state. In these methods the cluster pertur- 
bation theory is used that combines the exact diagonalization treatment of the multiband 
pd-model Hamiltonian inside the unit cell, and perturbation treatment of the intercell hop- 
ping in the X-operator representation. For the normal state this mean field approach is 
just a cluster generalization of the Hartree-Fock approximation for Hubbard fermions. For 
superconducting state, the mean field theory has been developed in the X-operator represen- 
tation that is reliable in the strong correlation regime 28j . Double occupation is prohibited 
in this approach by the local constraint formulated in the X-operator representation simi- 
lar to the local constraint in the slave boson approach. Contrary to the slave boson mean 
field theory where the local constraint is violated, in our mean field theory the X-operator 
algebra provides the local constraint in all stages of calculations. 

II. ISOTOPE EFFECT 

In the BCS-type approach to superconductivity a spin singlet pairing of the Hubbard 



fermions is given by the anomalous average ! 



-B q = (X^qX*' 2 ). For the d x 2_ y 2 -pairing 
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the gap equation reads [26] 

A k = ^ E J + \9 (le, - p| - Ud)}^ tanh (2) 

where r = /cbT, &b is Boltzmann constant, and T is temperature, y? q = (cosg^a — cosq y a)/2 
is the angle-dependent part of the gap A q = A </? q , the total coupling parameter in brackets 
is given by a sum of magnetic J and phonon A couplings. The ^-function as usually means 
that phonon pairing occurs in a narrow energy window of the ujd width near the Fermi energy. 
The normal phase dispersion £ q takes into account the spin correlation function c q , and 
three-cite interaction, the chemical potential \x is calculated self-consistently for the carrier 
concentration x in La2- x Sr x Cu04. The phonon coupling parameter A = f (x) G, where 
dimensionless function / (x) is given by f(x) = (1 + x) (3 + x)/8 — 3c i/4, and the parameter 
G is determined by the bare EPI matrix elements G = {gluck/^buck — 9 breath /^breath)- 

The concentration dependence of both magnetic and phonon couplings stems from the 
unusual statistics of the Hubbard fermions. Contrary to the free electron band with two 
electrons per atom the Hubbard subbands have the odd number of states which depends 
on concentration via total number of holes = 1 + x and the nearest neighbour spin 
correlation function cq\. Due to the antiferromagnetic type of correlation c i is negative 
and its contribution to A is positive. The appearance of the spin correlation function in 
the phonon coupling means some interference of the magnetic and phonon mechanisms of 
pairing. This function c 01 = 2 (SqS*) = ^SqSi^ characterizes the spin liquid properties of 
the underdoped cuprate and is concentration dependent [181 ] . 

All parameters in the gap equation §2§ but G have been obtained within ab initio 
LDA+GTB approach. The precision of modern calculations of the EPI matrix elements es- 
pecially for strongly correlated electrons seems to be not enough. Thus, different approaches 
give opposite conclusions. The larger EPI of the half-breathing mode versus buckling one has 



been discussed in papers 
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29 



30j. The average over Brillouin zone value A for buckling 



mode was obtained much larger than for breathing/half-breathing one 3l|. The ab initio 
calculations of the Eliashberg function a 2 F for YBa 2 Cu30y have revealed rather small EPI 



parameter A = 0.27 in disagreement with earlier approximate treatments 32|. Thus in this 



paper we consider the parameter G as the only fitting parameter, to find the value of G we 
calculate the isotope effect exponent determining as ao = — ^^j^ • Using the equation (J2J) 
we get the equation for the superconducting transition temperature T c then the analytical 
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expression for ao can be written in the form ao = where 
Intl = <f £ i*X ph 5 (|f q - fi\ - u D ) tanh 



(3) 

^2 = i E cosh" 2 {i^ J + A p/l (|£ q - fi\ - cu D )} 

The calculated oxygen isotope effect exponent ao as function of hole concentration is 
shown in the Fig. [TJ We have found that the positive (negative) sign of G results in positive 
(negative) sign of the exponent ao- The value G/J = 0.35 provides ao = 0.06 close to 
the La2- x Sr 3; Cu04 experimental data at the optimal doping. The increase of the isotope 
exponent away from the optimal doping is obtained. The higher value of the G parame- 
ter would result in the larger ao above the BCS value 0.5 but the optimal doping value 
ao (x = 0.17) will also increase. The doping dependence of the critical temperature T c is 
shown in the Fig. [2(a). It is clear that this dependence reproduces well the general structure 
of the superconducting dome with the optimal doping at x — 0.17 and disappearance of 
superconductivity in the underdoped region below x = 0.06. In the overdoped region the 
decrease of T c is more smooth than in experiment. Our approach from the undoped regime 
becomes less accurate in the overdoped region then at small doping. We have obtained that 
with positive parameter G the phonon contribution to pairing increases the magnetic one. 
The phonon contribution to the critical temperature is a little bit less than magnetic one but 
of the same order of magnitude. The absolute value of T c is too large and this is the general 
drawback of the mean field theory. The gap amplitude Ao as function of the hole concen- 
tration is shown in the Fig. Mjo)- Similar to the Fig. [2(a) the phonon contribution increases 
the gap value. At optimal doping the ratio 2Ao/kBT c = 4.9 is close to the experimental 
data 



We would like to discuss the effect of the strong EPI with apical oxygen breathing mode. 
This interaction is poor screened. The direct experimental proof of its importance is demon- 
strated by the colossal lattice expansion along c-axis in L^CuO^ under femtosecond inten- 



sive light pulses [34j . Nevertheless this large EPI does not contribute to the superconducting 
d-pairing due to the orthogonality of the in-plane electron momentum and c-axis phonon 
wave number for the apical oxygen breathing mode 26j |. The site selective isotope substi- 
tution confirms the absence of the isotope effect when the isotope is in the apical oxygen 

oxygen breathing mode contribution 



position 



35 



36(| . Previously the absence of the apica 



to the d x 2_ y 2 pairing has been obtained in the paper 
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FIG. 1: The doping dependence of the oxygen isotope exponent for the effective EPI parameter 
G/J = 0.35. 
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FIG. 2: The critical temperature (a) and the amplitude of the superconducting gap (b) as function 
of the doping concentration in La2-xSr x Cu04 for pure magnetic (G = 0) and both magnetic and 
phonon mechanisms (G/J = 0.35) of pairing. 

The effect of strong electron correlations on the phonon-mediated superconductivity and 
oxygen isotope effect has been considered in the slave-boson approach to the Hubbard 
model 37j] . The total agreement of our and slave-boson approaches is not possible. Thus 



the unphysical result of the zero bandwidth and zero EPI in the undoped limit 3?J is absent 
in our GTB method which provides the dispersion of the quasiparticles at the top of the 
valence band in undoped cuprates in agreement with the ARPES experiments. Nevertheless 
close to the optimal doping we can compare both approaches. Our EPI parameter increases 
with hole concentration x, and the X p d in the paper 37J also increases with x. The Coulomb 
interaction jj* SB our magnetic coupling which is the effect of strong 



correlations a 
in the paper 



so decreases with x. The isotope effect value has similar doping dependences 



37j ] and in our work. The phonon-mediated s-wave superconductivity with 



similar doping dependent isotope effect have been also discussed in the papers |38L I3QI] . 
The proper doping dependence of the oxygen isotope exponent a has been obtained 
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recently in the anharmonic version of the buckling mode mediated d x 2_ y 2 superconductiv- 
ity 40j. However this model neglects the depairing effect of the breathing and half-breathing 
mode and neglects the strong electron correlations which result in the magnetic mechanism 
of pairing. 

Summarizing our discussion we want to emphasize that both magnetic and phonon mech- 
anisms of d,^^ -pairing should be considered in realistic theory of superconductivity in 
cuprates. Many authors have previously discussed separately the magnetic mechanisms of 
pairing generated by strong electron correlations or phonon pairing which explains the iso- 
tope effect. Here we have shown that both mechanisms may work together increasing each 
other. Our theory is almost parameters- free. It is based on the correct microscopic descrip- 
tion of the undoped insulator state of La2Cu04 and the doping evolution of the emerging 
hole Fermi surface. All parameters of the electronic structure and the magnetic mechanisms 
of pairing have been calculated within the ab initio LDA+GTB approach. The only pa- 
rameter entering our theory is the combination of bare electron-phonon matrix elements G. 
Its sign G > is required to have the positive oxygen isotope exponent a, its value can be 
fitted to get the proper concentration dependence of a (x). The ab initio calculation of the 
electron-phonon matrix elements in the regime of strong electron correlations still remains 
the important unsolved problem. 
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